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Abstract 

Dynamical symmetry breaking in an expanding nuclear system is 
investigated in semi-classical and quantum framework by employing a 
collective transport model which is constructed to mimic the collective 
behavior of expanding systems. It is shown that the fluctuations in 
collective coordinates during the expansion are developed mainly by 
the enhancement of the initial fluctuations by the driving force, and 
that statistical and quantum fluctuations have similar consequences. 
It is pointed out that the quantal fluctuations may play an important 
role in the development of instabilities by reducing the time needed 
to break the symmetry, and the possible role of quantal fluctuations 
in spinodal decomposition of nuclei is discussed. 
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1 Introduction 



Dynamical symmetry breaking is a general behavior in nature and is observed 
in many field of physics. For example, phase transitions are often related to 
a dynamical symmetry breaking. The liquid-gas phase transition in infinite 
matter may be viewed as breaking of translational invariance which is related 
to mechanical instability of the system against small density fluctuations 
inside the spinodal zone. The spinodal instability may provide a possible 
mechanism for triggering the fragmentation of a hot piece of nuclear matter 
produced in heavy-ion collisions ffl, . The composite system formed in the 
overlap zone of the colliding ions expands and cools down. Depending on 
the initial compression and temperature, it may enter the unstable region of 
the phase diagram. In such a scenario, the fragment formation takes place 
by a rapid growth of density fluctuations in the spinodal region. Such a 
scenario is valid for infinite matter or a finite system. In matter, unstable 
modes are sound modes characterized by plane waves, whereas in a finite 
system these are collective vibrations associated with density fluctuations 
with multipolarity L and M. For example, a spherical system can break into 
pieces as soon as the collective mode with multipolarity L larger than one 
becomes unstable. 

Stochastic one-body transport models provides a suitable framework for 
a theoretical description of the collision process all the way from the initial 
stage up to the final state involving development of large density fluctua- 
tions and formation of fragments []3], ^, ||. This approach in a semi-classical 
approximation for the phase-space density, which is frequently referred as 
the Boltzmann-Langevin model, has been applied to investigate the spinodal 
decomposition of nuclear systems ||. Stochastic Time Dependent Hartree 
Fock (STDHF) approaches are now under development in order to take into 
account the quant al nature of single particle motion [|7|, §. 

A valuable insight on dynamical symmetry breaking, and the role of quan- 
tal and statistical fluctuations on the instabilities can be gained by investi- 
gating the early development of density fluctuations. The linear response 
treatment of the one-body stochastic transport models, in semi-classical or 
quantal framework, provides a useful basis for this purpose. In a recent work 
0, the dominant collective modes in an expanding finite nuclear system 
are identified in a linear response treatment by solving temperature depen- 
dent quantal RPA equation, and the properties of these collective modes 
are investigated as a function of dilution of the system. These calculations 
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illustrate that, as the system expands, the collective modes which were all 
stable at normal density, become softer until their energies become zero and 
then bifurcate toward the imaginary values, signaling the instability of the 
corresponding mode. 

The density fluctuations associated with the collective degrees of freedom 
can be determined in terms of collective transport equations. These transport 
equations for the collective variables can be deduced from the microscopic 
transport models by projecting the equation of motion onto a specific col- 
lective mode |nj [11], or can be derived by some other means ||12|| . Since, 
the standard one-body transport approach provides a classical description of 
the collective motion , the deduced collective models are best suited at high 
excitation energies when the temperature is much larger than the typical fre- 
quencies of the collective motion. On the other hand, when the characteristic 
energies of the collective modes are comparable with or larger than tempera- 
ture of the system, the quantal fluctuations of the collective dynamics become 
important or even dominant. Then, a quantum transport approach of the 
collective variables is required for describing dynamics of large density fluc- 
tuations. The stochastic one-body transport description can be improved by 
incorporating the memory effects associated with finite duration of binary 



collisions. As demonstrated in a recent work ||13| , |14}| , in such non-Markovian 
description, transport properties of collective motion evolve in accordance 
with the quantal fluctuation-dissipation relation |L5] . 



In this work, we consider a single collective mode, specifically the quadru- 
pole moment of the fluctuating density in an expanding nuclear system and 
study its evolution by employing collective transport models, in both clas- 
sical and quantal forms. Although these models can be derived from the 
underlying microscopic transport models in Markovian and non-Markovian 
forms, we follow a phenomenological approach and parameterize the basic 
ingredients of the model in accordance with the recent TDHF calculations 
for expansion of a hot spherical system |jl6|, O, and the RPA calculations 



for describing the unstable collective modes 0. In the Fokker-Planck ap- 
proach, the coupling between the collective mode and the intrinsic degrees 
of freedom is described in terms of friction and diffusion coefficients. This 
simple but realistic model allows us to calculate the development of fluctua- 
tions associated with the collective mode during the expansion of the nuclear 
system, and to investigate the role of quantal and statistical fluctuations in 
the dynamics of the symmetry breaking. 
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2 The model 



2.1 The collective Hamiltonian 

Combining the results of the finite temperature TDHF description of the 



expansion of hot nuclear sources ref. fig , 17 and the RPA investigation of 
the stability of a finite hot expanding nuclear source ||, we can determine 
the most important collective modes and describe their transitions from sta- 
ble to unstable configurations as a function of expansion, temperature and 
source size. Here, we specifically consider evolution of the quadrupole shape 
fluctuations in an expanding spherical hot nuclear source. According to fi- 
nite temperature TDHF calculations with initial temperatures in the range 
of T = 5 — 10 MeV, the system executes monopole vibrations and cools down 
by particle evaporation. Depending on the size of the source and the initial 
temperature, it takes between 30 fm/c and 75 fm/c to reach the maximum 
dilution that occurs typically around 1/3 of the normal nuclear matter den- 
sity. According to the constraint RPA calculations, the quadrupole mode 
with an initial frequency cu, which is found to be about u ~ QhA^ 1 /^ MeV % 
at moderate temperature, becomes softer as the system expands. It eventu- 
ally becomes unstable and reaches an imaginary frequency of about i20A~ 1 / 3 
MeV/h at the maximum dilution. 

In order to describe the evolution of the quadrupole mode during the ex- 
pansion of the source, we consider the following model, in which the Hamil- 
tonian of the collective motion is taken as 

H=^ + U(Q) (1) 

where M = 3AmRo/87r is the irrotational mass parameter associated with 
the collective variable Q determined by the sharp nuclear surface Rq, A is 
the number of nucleons in the source and U(Q) represents the collective 
potential energy. For small deformation Q , the potential energy has a har- 
monic form U (Q) = K (t) Q 2 /2 + ... with a time dependent spring constant 
K (t) = K g(t) where K = uo 2 M and uo = MeV/h denotes the 

frequency of the giant quadrupole vibrations at normal density. The time 
dependent factor g(t) in the spring constant describes the transition from 
stable to unstable regime during the expansion of the source. This factor is 
parameterized In accordance with the TDHF and RPA calculation: starting 
from g(0) = 1 it monotonicaly decreases in time and reaches a value around 
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9{to) = —1/9 at a characteristic time to f° r the source to reach the maximum 
instability (typically, t = 50 fm/c). The harmonic approximation is valid 
during early stages of the unstable evolution when the magnitude of fluctua- 
tions is small. In order to describe large fluctuations, we need to incorporate 
an anharmonic terms in the potential energy, U (Q), which prevent the fluc- 
tuation to grow indefinitely without bound. For symmetry considerations, 
we add an anharmonic term of the form Q 4 in the potential energy to achieve 
the saturation of the fluctuations 

U(Q) = ^Q 2 + PQ* (2) 

Since, we are mainly concern with the early development of fluctuations, the 
precise form of this non-linear term is not important for this purpose. The 
numerical coefficient in front of the non-linear term is taken to be f3 = K /2 
in order to obtain saturation of fluctuations at a reasonable value of the 
collective variable Q which corresponds to the quadrupole deformation of 
the fragmenting system (e.g. the quadrupole deformation may lead to break- 
up of a system into two pieces). The numerical calculations presented are 
performed for a system containing A = 64 nucleons. 



2.2 Classical Langevin dynamics 

The classical Langevin approach provides a useful framework for describing 
the dynamics of the collective motion including dissipation and fluctuations. 
In this approach, the trajectory in the collective phase space, (g»(t),p«(t)) 
of an event i, is determined by solving a Langevin equation, 



dt M 
d pii) dU (Q« 



dt dQ 



- 7 P« + £«> (t) (4) 



where 7 is the friction coefficient and £^(t) is the random force arising from 
the coupling of the collective degrees of freedom with the single particle 
motion. The fluctuating force £®(t) is assumed to be a white noise, and 
specified by a Gaussian distribution with zero mean < £(£) >= and a 
second moment 

(Z(t)ZV)) = 8(t-lf)2D (5) 
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where (...) represents the statistical average over the ensemble of events i. 
The diffusion coefficient D is determined in terms of the friction coefficient 
by the classical dissipation-fluctuation relation 

D = jMT (6) 

where T represents the time-dependent temperature of the internal system 
which may be specified by assuming an isentropic expansion of the nuclear 
source. The only remaining parameter of the model is the friction coeffi- 
cient, which may be estimated from the observed damping width of the giant 
quadrupole excitation according to 7 = T/%. In the calculations, we take the 
damping width to be T = 85A~ 2 / 3 MeV. 

By generating a sufficiently large number of events, we can construct the 
phase-space density f(P,Q,t) associated with the collective variables P and 
Q as 

1 N 

f(Q,P,t) = -J2$(Q- Q (l) (t)) 5 (P ~ P {l) (t)) (7) 

8=1 

where N denotes the number of events. The probability distribution n(Q,t) 
of Q and the variances of the collective variables ctq(£) and crp(t) are calcu- 
lated as 



r 1 N 

n(Q,t) = J dP f(Q,P,t) = -Y,6(Q-QM(tj) (8) 



and 



a Q (t) = <Q 2 >= J dPdQQ 2 f(Q,P,t) (9) 
CTp(t) = < P 2 >= J dPdQ P 2 f(Q 1 P,t) (10) 

The initial conditions of the events generated by the Langevin equation is 
specified by the Boltzmann distribution with the initial temperature T 

/b(M = I-p(-M). (id 

In our simulations, we sample the initial fluctuations using a metropolis 



algorithm[lS 



The classical description of the Langevin equation provides a good ap- 
proximation at sufficiently high temperatures at which the dynamics is dom- 
inated by thermal fluctuations. At low temperatures, the Langevin approach 
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may be improved by incorporating the quantal fluctuations through the ini- 
tial conditions and by specifying the diffusion coefficient through the quan- 



tal fluctuation-dissipation relation |13|, OJ]. As a matter of fact, this de- 
scription produces exact quantal phase-space density in the harmonic limit, 
for which the quantal fluctuation-dissipation relation may be expressed as 
D = ( , -fMh,u)/2)coth.(hij/2T) for real u. However, in general for non-linear 
evolution, collective motion should be investigated in the basis of a quantal 
transport equation. 

2.3 Quantal transport model 

In quantal approaches the system is described by its density matrix pit) on 
the collective Hilbert space, p(Q,Q',t) = (Q\ p{t) \Q') , which is determined 
by a quantal transport equation 

ifi^p(t) = [H(t),p(t)] +ihK(p) (12) 

where K(p) represents a "collision term" arising from the coupling of the 
collective motion with the intrinsic variables, which corresponds to the fric- 
tion and diffusion term in the Langevin equation. The collective density 
matrix may be expanded in terms of a complete set of time-dependent wave 
functions \^i(t) > as 

p(t) = j2\Mt))Pi j (t)(s> j (t)\ (is) 

where the time-dependent wave functions are obtained by solving the time- 
dependent Shrodinger equation 

ih^ \%(t)) = H(t) \%(t)) (14) 
and the elements of density matrix pij (t) obey a master equation 

^p ij (t)=<%(t)\K(p)\V j (t)> (15) 

If the system initially is in statistical equilibrium at an initial temperature 
T , it can be represented by a density matrix, 

^ = °)= <h x p(-|0 (16) 
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It is convenient to introduce the eigenstates \ipi > and the eigen-energies 
of the initial Hamiltonian H = H(t = 0) 

H \(pi > = Ei\(pi > (17) 

as initial basis. Consequently, the initial density matrix reduces to Boltz- 
mann occupation probabilities of the eigenstates of the initial Hamiltonian, 

1 / E\ 
Pij(t)=6 ij -exp --f . (18) 



^ V T 

The probability density n(Q,t), and the variances cq(£) and ap(t) of the 
collective variables are calculated according to 

n(Q,t) = {Q\p(t)\Q) (19) 

and 

cr Q (t) =< Q 2 >= tiQ 2 p(t) (20) 

a P (t) =< P 2 >= trP 2 P (t) (21) 

Since the model Hamiltonian is nearly harmonic at the initial time t — 0, the 
initial variances of the collective variables can approximately be expressed in 
a closed form as 

o-q(O) = (Q 2 )o - (hcu/2K Q )coth(huj/2T ) -> hw/2K Q 

Io ^° (22) 
^(0) = (P 2 ) ~ (Mhu/2)coth(hu/2T ) -» M^/2 

Here, the approximate expressions are valid when temperature is small as 
compared to the collectives frequency, To fca;, and hence the initial vari- 
ances are dominated by the ground state zero point fluctuations. In the 
opposite limit when temperature is much larger than half the initial collec- 
tive energy htu 2T the initial variances of the collective variables can be 
approximated by 

a Q (0) ~ T /K 

a P (0) ~ MT 1 6) 

which corresponds to the classical description discussed in the previous sec- 
tion. 



8 



3 Results 



We carry out a number of simulations on the basis of two different approaches 
presented above to investigate the development of quadrupole shape instabil- 
ities in an expanding spherical nuclear system prepared at a range of initial 
temperatures and calculate the probability density n(Q,t) and the variances 
<7q(£), <Jp{i) of the collective variables. In our calculations, we consider that 
the initial state is prepared at normal density and at a temperature in the 
range of T = 5 — 10 MeV. This is in accordance with the BUU calcula- 
tions for central heavy-ion collisions around Fermi energy, which show that 
the collisions lead to an equilibrated composite spherical system at tempera- 
tures around T = 10 MeV. The second important parameter in our model is 
the characteristic time to for expansion until the system reaches the turning 
point, which has a typical magnitude around to = 50 fm/c. In order to see 
the influence of the expansion time on the dynamics of symmetry breaking, 
we present calculations for different t . 

3.1 Classical dynamics 

The result of simulations of the classical Langevin equation is illustrated in 
figure 1 by solid lines, in which the probability distribution n(Q,t) and the 
potential energy U(Q,t) are plotted as a function of Q at several times for 
two different initial temperatures T = 5 MeV (left panel) and T = 10 MeV 
(right panel). In these simulations temperature is assumed to be constant 
and equal to the initial value. The width of distribution rapidly increase 
as the potential gets softer and the probability distribution splits into two 
components after the potential exhibits a local maximum at Q — 0. Figure 2 
displays the variances associated with the collective variable Q and the col- 
lective momentum P as a function of time for the initial temperatures T = 5 
MeV and T = 10 MeV. The fluctuation in the collective coordinate rapidly 
increase until the system reaches the maximum instability and then saturates 
after 100 fm/c. During the same time the momentum distribution remains 
almost constant and goes to the asymptotic value given by the equipartition 
theorem o P ~ MT . In figures 1 and 2, the dotted lines display the results 
of the pure potential calculations without the dissipation and fluctuations. 
During the expansion phase until the source reaches maximum dilution at 
about t = 50 — 60 fm/c, the pure potential result for n(Q,t) is very close to 
the Langevin simulations with friction term and stochastic force. Therefore, 
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Figure 1: Probability distribution of collective variable and potential energy 
as a function of Q at times t = 0, 50, 100 fm/c for initial temperatures To = 5 
Me V (left panel) and To = 10 MeV (right panel). Solid, dashed and dotted 
lines are the Langevin simulations with constant initial temperature, the 
Langevin simulation with time dependent temperature and the pure potential 
calculations, respectively. 
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the fluctuations of the collective variable Q during the expansion phase are 
not developed dynamically but originate mainly from the propagation of the 
initial statistical fluctuations by the driving force. As the system expands, 
it cools down and consequently the magnitude of the stochastic force in the 
Langevin equation become smaller. The dashed lines in figures 1 and 2 indi- 
cate the simulations carried out by employing a time dependent temperature, 
which is determined according to T (t) = To(R /R(t)) 2 by assuming adia- 
batic expansion, where R (t) represents the root-mean-square radius of the 
expanding system. As seen from these, the cooling has a minor effect on 
the probability distribution n(Q, t) during the early phase of the expansion. 
However, cooling has a sizable effect in momentum space at large times, as 
seen in the asymptotic value of the op which also come close to the result 
obtained by the pure potential calculations. 

3.2 Comparison between quantum and classical dy- 
namics 

According to the classical Langevin simulations, during the expansion phase, 
the fluctuation of collective variable Q is mainly determined by the potential 
evolution and the friction and stochastic forces have a minor effect. It is 
reasonable to expect a similar behavior in the quantal evolution of the col- 
lective motion. Therefore in the application of the quantal model presented 
in section 2.3, we neglect the collision term, and calculate the probability 
distribution n(Q,t) and the variances <?Q{t), &p{t) by keeping the occupa- 
tion probabilities pi to be constant and equal to initial values specified by 
the Bolztmann factors. Solid lines in figure 3 illustrates the result of quan- 
tal calculations of the probability distribution n(Q,t) at different times for 
two initial temperatures T = 5 MeV (left panel) and T = 10 MeV (right 
panel). In the same figure, dotted lines show the classical simulations with- 
out friction and stochastic force, and also, as a reference dashed lines show 
the quantal calculations performed only with the ground state (ie the T = 
MeV case). Figure 4 illustrates the evolution of the variances &Q(t) and 
<jp{t) as a function of time for initial temperatures T = 5 MeV and T = 10 
MeV. The classical simulations, the full quantal calculations and the quantal 
calculations with ground state are shown by dotted, solid and dashed lines, 
respectively. The fluctuations of the collective variable, 0"q(£), exhibit an 
exponential growth as soon as the mode becomes unstable, and saturate at 
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Time (fm/c) 



Figure 2: Variances associated with collective variable Q and collective mo- 
mentum P as a function of time at initial temperatures T = 5 MeV(left 
panel) and T = 10 MeV (right panel). Solid, dashed and dotted lines are 
the Langevin simulations with constant initial temperature, the Langevin 
simulation with time dependent temperature and the pure potential calcula- 
tions, respectively. 
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Figure 3: Probability distribution of collective variable and potential energy 
as a function of Q at times t = 0, 50, 100 fm/c for initial temperatures To = 5 
Me V (left panel) and To = 10 MeV (right panel). Solid, dashed and dotted 
lines are the full quantal calculation, the quantal calculations with ground 
state, and the classical calculations without friction and stochastic forces, 
respectively. 
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Figure 4: Variances associated with collective variable Q and collective mo- 
mentum P as a function of time at initial temperatures To = 5 MeV(left 
panel) and T = 10 MeV (right panel). Solid, dashed and dotted lines are 
the full quantal calculation, the quantal calculations with ground state, and 
the classical calculations without friction and stochastic forces, respectively. 
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large times, in both the classical and quantal calculations. The exponential 
growth rates are the same in both calculations, but the magnitude of fluc- 
tuations are larger in the quantal case. This is illustrated in Figure 5 which 
shows the critical time r it takes to reach oq = 0.2 as a function of the 
expansion time to of the source for the initial temperatures T = 5 MeV and 
T = 10 MeV. As seen, the quantal fluctuations shown by dashed lines reduce 
the time it takes to reach a finite value of the fluctuations. The reduction 
of this time, depending on the temperature of the system, can be as large as 
factor of two. For comparison, figure 6 shows different classical calculations 
of the critical time r as a function of the expansion time to. In this figure, the 
Langevin simulations with constant initial temperature, the Langevin simu- 
lations with time dependent temperature and the pure potential calculations 
are indicated by solid, dashed and dotted lines, respectively. 

4 Conclusions 

We investigate the development of density fluctuations associated with collec- 
tive modes (here we specifically consider quadrupole mode) in an expanding 
nuclear system by employing phenomenological classical and quantal trans- 
port equations. According to mean-field calculations, a finite nuclear system 
prepared at moderate temperatures around T = 5— 10 MeV expands until the 
turning point which is located inside the unstable zone, and then continue to 
execute monopole vibrations. We calculate the evolution of the probability 
distribution of the collective coordinate Q and the collective momentum P 
and determine the variances associated with these distributions as a function 
of time. Our investigations indicate that during the expansion phase, the dis- 
sipation and fluctuation mechanism play an important role in the evolution 
of the momentum distribution, but the fluctuations in collective coordinate Q 
are developed mainly by enhancement of the initial fluctuations of statistical 
or quantal origin. Furthermore, it appears that quantal fluctuations, which 
are already present at the initial state, play an important role in particular at 
low temperatures, in dynamics of symmetry breaking by reducing the time 
required for reaching the critical fluctuations. Consequently, the spinodal de- 
composition and the resultant fragmentation may take place in a faster time 
scale than predicted by semi-classical simulations || Therefore, it is off 
great interest to develop stochastic simulation methods for heavy-ion colli- 
sions by incorporating quantal fluctuations associated with collective motion 
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Figure 5: Time r that takes to reach oq = 0.2 as a function of the expansion 
time to at temperatures T = 5 MeV (top panel) and T = 10 MeV (bottom 
panel). Solid, dashed and dotted lines are the full quantal calculation, the 
quantal calculations with ground state, and the classical calculations without 
friction and stochastic forces, respectively. 
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Figure 6: Same as in figure 5. Solid, dashed and dotted lines are the Langevin 
simulations with constant initial temperature, the Langevin simulation with 
time dependent temperature and the pure potential calculations, respectively. 
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in a suitable manner. 
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